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We make a Monte Carlo study on compositeness of first generation quarks and leptons using the 
Drell-Yan distribution in the high dielectron mass region at the Tevatron and LHC energies. The 
current experimental lower limits on the compositeness scale, A, vary from 2.5 to 6.1 TeV. In the 
present analysis, we assume that there will be no deviation of the dielectron mass spectrum from 
Standard Model prediction at center of mass energy 2 TeV (Tevatron) and 14 TeV (LHC). We then 
find that in the LL, RR, RL and LR chirality channels of the quark-electron currents, it is possible 
to extend the lower limits on A (at 95% CL) to a range of 6 to 10 TeV for 2 fb' 1 and 9 to 19 TeV 
for 30 /ft -1 of integrated luminosity at Tevatron. At LHC, the corresponding limits extend to a 
range of 16 to 25 TeV for 10 fb' 1 and 20 to 36 TeV for 100 /6" 1 of integrated luminosity. 

PACS numbers: 12.60.Rc, 12.60.-i, 13.85.-t 
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The proliferation of quarks and leptons has inspired 
the speculation that they could be composite structures, 
i.e. bound states of more fundamental constituents of- 
ten called preons . Below a characteristic energy scale 
called the compositeness scale, A, the preon-binding in- 
teraction becomes strong and binds the constituents to 
form composite states like the quarks and leptons. With 
such a composite structure, there would be significant 
deviation from the Standard Model (SM) prediction of 
high energy cross sections. No such deviation has been 
observed so far. These null results have been used to 
put lower limits on quark-lepton compositeness scale A. 
which varies from 2.5 to 6.1 TeV || in the various chi- 
rality channels of the quark-lepton currents. 

In this paper, we consider the effects of composite 
structure of first generation quarks and leptons on the 
Drell-Yan (DY) process qq — > e + e~ HJ. If the composite- 
ness scale, A, is much greater than Vs, the center of mass 
energy of the colliding partons, the quarks and electrons 
would appear to be point-like. The substructure coupling 
can then be approximated by a four-fermion contact in- 
teraction giving rise to the following effective lagrangian[] 



g 2 

Cql = f?LLfe7 M 9z,)(eL7M e i) + r ?£-R(9L7 M 9L)(eii7M e -R) (0-1) 

where 



is the left-handed quark doublet, ur and da are the 
right-handed quark singlets, ez, and en are the left- 
and right-handed electrons respectively. The composite- 
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ness scale (A) is chosen so that %l = 1 and the largest 
| rjij |= 1, where go is the coupling constant for the con- 
tact interaction and rjij is the interference term between 
the contact interaction and the SM lagrangian for the 
ij th channel, with i and j representing the helicities of 
the quark and the lepton currents. Including the above 
contact interaction (at A >> \/§), the DY cross section 
gets transformed as H: 



dm dm 



(0.2) 



where j3 — 1/A 2 and m is the dielectron invariant mass. 
In this expression, / is due to the interference of DY and 



1 Here we have assumed that the contact interaction is color singlet 
and weak-isoscalar. 



the contact term, and C is the pure contact term contri- 
bution to the cross-section. The deviation in the dielec- 
tron production from SM expectations would be domi- 
nant in the high mass region above the Z pole. We have 
made separate studies for quark- electron compositeness 
for an integrated luminosity of 2 fb^ 1 (Run II) and 30 
fb^ 1 (TEV33) with respect to the D0 detector at Teva- 
tron and an integrated luminosity of 10 fb^ 1 and 100 
fb^ 1 with respect to the CMS detector at LHC. How- 
ever the results should be valid for the CDF detector at 
Tevatron and the ATLAS detector at LHC as well. Wc 
have simulated dielectron production through DY pro- 
cess alone in pp (pp) collisions at center of mass energy, 
y/s, equal to 2 TeV (14 TeV) using PYTHIA §. How- 
ever since PYTHIA does not incorporate all the com- 
positeness models, we have used a separate parton level 
Monte Carlo program to estimate dielectron production 
rates in the presence of compositeness. Assuming that 
the Tevatron and LHC data on dielectron production are 
consistent with DY predictions under SM, we extract lim- 
its on compositeness scale using Bayesian technique of 
statistical inference We have considered four dif- 

ferent models corresponding to the LL, RR, RL and LR 
chirality channels of equation 0.1 for quark-electron com- 



positeness. The choice of rjij for the different models of 
compositeness is listed in Table ffl. 



TABLE I. Choice of ijij for different contact interaction 
models. The superscript on the model denotes the nature 
of interference between the contact interaction and the SM 
lagrangian. Constructive interference (77^ = — 1) is denoted 
by a + and destructive interference (rjij = + 1) is denoted 
by a -. 



Model 


f]LL 


Vrr 


r/LR 


r\RL 


LL± 


Tl 











RR ± 





Tl 








LR ± 








Tl 





RL± 











Tl 



Exploring the lower limits on A at Tevatron 

We simulate pp collisions using PYTHIA at 2 TeV and 
generate DY dielectron events between 95 GeV and 1.5 
TeV of the dielectron invariant mass. The total number 
of dielectron events generated by PYTHIA, N gen , gives 
the expected number of background subtracted dielec- 
tron events, Nuy, to be collected at Tevatron as : 



N, 



DY 



= e x N, 



gen 



(0.3) 



where e is the detection efficiency of the dielectron. The 
detection efficiency, e, of the dielectron involves contri- 
bution from the following terms: 

(a) Energy smearing, 
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(b) Electron identification efficiency, ei, and 

(c) Acceptance, t%. 

The energy resolution of the electromagnetic calorimeter 
of the upgraded D0 detector is parameterized as : 



(J) 2 = + ( 



f (E in GeV) 



(0.4) 



where the constant term, C, and the stochastic term, a 
are taken to 2 % and f6 % respectively . We take the 
electron identification efficiency, e%, for a single electron 
to be 85%. The identification efficiency for a dielectron 
is then e\. The acceptance, £2, of dielectron events in pp 
collisions is defined as the fraction of events in which the 
e + e~ pair passes the fiducial and the kinematic cuts after 
taking into account the energy smearing. The fiducial 
and the kinematic cuts used are: 

• I T) |< 2.5, where r\ is the pseudorapidity (= 
— ln[tan(^)]). This ensures that the dielectron 
event selected is in the active detector region. 

• A kinematic cut of px > 25 GeV, where pt is 
the transverse momentum of the electron and the 
positron. This cut ensures an efficient trigger^] 

The dielectron detection efficiency, e, is then : 



(0.5) 



We then generate the expected number of dielectron 
events, N^ xp , in various mass bins including the effect 
of the composite structure of quarks and electrons for 
various values of A using the parton level Monte Carlo. 
We calculate the cross section (cr A ) for the production of 
dielectrons including t erm s from the contact interaction 
lagrangian of equation 3.1 with the SM lagrangian. The 
LO cross section calculation is corrected for higher order 
QCD effects using a K-f actor of 1.220. We checked the 
parton level MC calculation by comparing its prediction 
with that from PYTHIA for the Drell-Yan process. Both 
calculations agree to within a few percent as shown in 
Fig. § 

In order to obtain the lower limit on A, we then use 
the Bayesian technique to compare the Drell Yan dielec- 
tron mass distribution (ie., Ndy) in the high mass re- 
gion with the expected dielectron mass distribution for 
various values of A (ie., N^ xp ). Limits are obtained inde- 
pendently for each separate channel of the contact inter- 
action lagrangian: LL, RR, RL and LR with rjij = ±1 . 
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FIG. 1. Dielectron invariant mass spectra between 80 GeV 
and 1.5 TeV for DY process at yfs = 2 TeV , as predicted 
by PYTHIA and as calculated using our parton level Monte 
Carlo. 



Fig. g shows the cross section versus the dielectron in- 
variant mass, in the high mass region between 50 GeV 
and 1.8 TeV in the LL channel for different values of A 
for rjij = — 1 (constructive interference) and Fig. || 
shows the corresponding plot for rjij = + 1 (destruc- 
tive interference). 
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FIG. 2. Cross section, Act (in 50 GeV bins), versus dielec- 
tron invariant mass, m, between 50 GeV and 1.8 TeV for DY 
process and three different values of A in the LL channel for 
rjij = - 1. 



2 This cut is based on the D0 Run I analysis of DY data at 1.8 
TeV §. 

3 This K-factor is the ratio of the NNLO DY cross section to the 
LO DY cross section at 1.8 TeV ^,|uj. We consider the same value 
for the K-factor for DY + compositeness at 2 TeV. 



Since the effect of compositeness is most pronounced 
in the high dielectron mass region we consider 10 differ- 
ent mass bins of variable width between 120 GeV and 1.5 
TeV. The expected number of events at the composite- 
ness scale, A, in the k th mass bin is given as : 
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FIG. 3. Cross section, Act (in 50 GeV bins), versus dielec- 
tron invariant mass, m, between 50 GeV and 1.8 TeV for DY 
process and three different values of A in the LL channel for 
VH = + 1. 



N. 



A,k 



e fc (a A ' k x L), 



(0.6) 



where a A ' k is the cross section (including compositeness) 
for the k th mass bin and L is the integrated luminosity. 
The posterior probability for the compositeness scale to 
be A given the expected DY dielectron data distribution, 
d , is: 

n 

P(A \do) = -= I] pk ( N DY I N e % k ) P(e k , L, A) (0.7) 



k = 1 

where Z is the normalization constant. P k (N^ Y I ^exp) 
is the likelihood function which follows a Poisson distri- 



bution for small N, 



P k (N k Y | Nt$) 



A,k. 
exp " 



n dy'- 



(N^ k < 10) (0.8) 



and a Gaussian distribution for large N^ k , 
N A ji k and standard deviation, o~\ , (ci = 



with mean 

A,k 
exp 
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P(AT| 



I N A,k ) 

I 1 " exp ) 



{Ntl > 10). (0.9) 



P(e k , L, A) is the joint prior probability for the dielec- 
tron detection efficiency, e k , the integrated luminosity, L, 
and the compositeness scale, A. Taking e k , L and A to 
be independent, 



P(e k ,L,A) = P{e k ) P{L) P{A). 



(0.10) 



The prior probabilities of detection efficiency, e k , and in- 
tegrated luminosity, L, are assumed to be Gaussian with 
their estimated value in each bin as the mean and corre- 
sponding error as the width of the Gaussian. The prior 
distribution P(A) is chosen to be uniform in 1/A 2 . This 



represents a prior essentially flat in cross section. The re- 
sulting posterior density P(A | do) peaks at 1/A 2 = 
and falls off monotonically with increasing 1/A 2 . The 
95% CL lower limit on A is defined by: 



dA' P(A' | d ) = 0.95. 



(0.11) 



The values of efficiency, e k , and the expected number 
of DY events, A^ y [| in individual mass bins arc listed 
in Table [fi] for an integrated luminosity of 2 fb^ 1 and 30 
fb~ l . The expected 95% CL lower limits on A for the LL, 
RR, RL and LR helicity channels of the quark- electron 
currents for both constructive and destructive interfer- 
ence are listed in Table II] and Table |v] for integrated 
luminosities of 2 fb^ 1 



and 30 fb respectively. 



TABLE II. Detection efficiency and expected number of 
DY events in different mass bins 



mass bin 


N k y 


N k Y 


(GeV) 


L = 2 fb- 1 


L = 30 fb- 1 


120-160 


0.590 2335.8 


34508.1 


160-200 


0.629 606.9 


8990.1 


200-240 


0.655 236.3 


3589.4 


240-290 


0.663 117.8 


1942.8 


290-340 


0.675 66.5 


877.1 


340-400 


0.668 34.0 


461.7 


400-500 


0.689 23.8 


276.0 


500-600 


0.712 6.5 


98.3 


600-1000 


0.677 1.5 


42.6 


1000-1500 


0.723 


2.2 


TABLE III. 


Expected 95% CL lower 


limits, Au m , on 


the compositeness scale for different helicity channels of the 


quark-electron currents for L = 2 fb 


1 at 2 TeV with 


5e k = 15 % and SL = 5 % 






Aum (TeV) 


A lim (TeV) 


Channel 


iVH = - 1) 


(Vij = + 1) 


LL 


10.1 


8.0 


RR 


9.3 


6.0 


RL 


7.8 


5.7 


LR 


7.3 


6.0 



Exploring the lower limits on A at LHC 

We have made a similar analysis of the DY process 
including the effect of quark-electron compositeness at 
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TABLE IV. Expected 95% CL lower limits, A lim , on 
the compositeness scale for different helicity channels of the 
quark-electron currents for L — 30 at 2 TeV with 



15 % and SL 



5 % 



Channel 


Au m (TeV) 

(vn = - 1) 


A iim (TeV) 

(m = + i) 


LL 


18.9 


17.8 


RR 


17.0 


15.1 


RL 


13.5 


9.1 


LR 


12.1 


9.2 



14 TeV. As before we have assumed that DY dielectron 
data that would be collected by the CMS detector at 
LHC would agree with SM prediction. We then use the 
Bayesian technique to obtain the lower limits on A at 14 
TeV. We have made separate studies for 10 fb- 1 of data 
and 100 fb' 1 of data. A K-factor of 1.13 §,[1 has been 
used as the NNLO correction factor. Fig. ishows the 
cross section versus the dielectron invariant mass, in the 
high mass region between 50 GeV and 2 TeV in the LL 
channel for different values of A for rftj = — 1 (construc- 
tive interference) and Fig. || shows the corresponding plot 
for rjij = + 1 (destructive interference). 

We generated DY events in the dielectron mass range 
of 150 GeV to 2 TeV. We then compared the expected 
number of DY events, Ndy, at y/s = 14 TeV with the 
expected number of dielectron events, , at various 
values of A in the mass range of 500 GeV to 2 TeV where 
the deviation from SM predictions due to the composite 
structure of quarks and electrons is most pronounced at 
LHC. The electron identification efficiency, e\, is taken 
to be 95% (jT2|. The constant and stochastic terms in 
the energy resolution of the electromagnetic calorimeter 
of the CMS detector are taken to be [fl2f : 



C = 0.55%, and 
a = 2.7%, | 77 | < 1.5 
5.7%, 1.5 < \rj 



(0.12) 



< 2.5 



The fiducial and kinematic cuts selected are the same 
as for D0. The values of e k and Y£, y [| in individual 
mass bins are listed in Table M for integrated luminosi- 
ties of 10 fb- 1 and 100 fb~\ The expected 95% CL 
lower limits on A for the LL, RR, RL and LR helicity 
channels of quark-electron currents for both constructive 
and d estr uctive interference are listed in Table VI and 
Table VII for integrated luminosities of 10 and 100 
fb~ x respectively. 

The discovery limits for A (defined as a deviation of 5a 
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FIG. 4. Cross section, Act (in 50 GeV bins), versus dielec- 
tron invariant mass, m, between 50 GeV and 2 TeV for DY 
process and three different values of A in the LL channel for 
rjij = - 1. 

from SM prediction) for the various models have been 
listed for integrated luminosities of 10 fb -1 , 50 /b " 1 , 
100 fb- 1 , 200 fb- 1 and 500 fb- 1 in Table |vm] for 
rjij = — 1 and in Table IX for 77^ = 



1. 



is generated with a K-factor of 1.13 in PYTHIA. 



Plots of the discovery limit versus the integrated lu- 
minosity for the various chirality channels are shown in 
Fig. H for rjij = — 1 and in Fig. ^ for rjij = +1. 

To conclude, we have performed a Monte Carlo study 
of the dielectron invariant mass spectrum (DY + com- 
positeness) for pp collisions at 2 TeV and pp collisions 
at 14 TeV. We have considered the LL, RR, RL and LR 
chirality channels of the quark-electron currents. Assum- 
ing that Standard Model will describe the high mass DY 
dielectron data at 2 TeV and 14 TeV we have found that 
it is possible to extend the lower limits on the compos- 
iteness scale, A, from the existing limits. 



• For pp collisions at Tevatron we have made sepa- 
rate studies for integrated luminosities of 2 fb- 1 
and 30 fb -1 with respect to the D0 detector. The 
expected 95 % CL lower limits on A range between 
6 to 10 TeV and 9 to 19 TeV for 2 and 30 fb' 1 
of dielectron data, respectively. These limits are in 
agreement with similar limits on A quoted between 
6 to 10 TeV for 2 fb- 1 and 14 to 20 TeV for 30 fb- 1 
of data with respect to the CDF detector at Teva- 
tron Q. 

• For pp collisions at LHC we have considered 
10 fb -1 and 100 fb -1 of dielectron data with re- 
spect to the CMS detector. The expected 95 % CL 
lower limits on A range between 16 to 25 TeV for 
10 fb' 1 and between 20 to 36 TeV for 100 fb- 1 of 
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FIG. 5. Cross section, Ait (in 50 GeV bins), versus dielec- 
tron invariant mass, m, between 50 GeV and 2 TeV for DY 
process and three different values of A in the LL channel for 
Vij = + 1- 

dielectron data. 



TABLE V. Detection efficiency and expected number of 
DY events in different mass bins 



mass bin 
(GeV) 


k 

e 


Ndy 
L = 10 fb^ 1 


Ndy 
L = 100 fb 


500-510 


0.660 


56.86 


586.63 


510-520 


0.617 


44.22 


506.30 


520-530 


0.668 


48.74 


493.67 


530-540 


0.654 


44.22 


425.98 


540-550 


0.647 


37.91 


418.76 


550-560 


0.662 


46.03 


378.15 


560-570 


0.666 


37.91 


361.90 


570-580 


0.668 


34.30 


329.41 


580-600 


0.673 


45.13 


552.33 


600-625 


0.684 


57.76 


621.82 


625-650 


0.678 


50.54 


509.91 


650-675 


0.681 


45.13 


419.66 


675-700 


0.692 


28.88 


386.27 


700-750 


0.717 


52.35 


589.33 


750-800 


0.728 


36.10 


452.15 


800-900 


0.731 


53.25 


613.70 


900-1000 


0.756 


31.59 


336.63 


1000-1200 


0.752 


36.10 


315.88 


1200-1400 


0.782 


16.25 


151.62 


1400-2000 


0.791 


16.25 


135.38 



• We have also explored the discovery potential for 
quark-electron compositencss (defined as a devia- 
tion of 5cr from SM prediction) at LHC as a func- 
tion of integrated luminosity. 
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TABLE VI. Expected 95% CL lower limits, A Um , on 
the compositeness scale for different helicity channels of the 
quark-electron currents for L = 10 /& 1 at 14 TeV with 
Se k = 15 % and SL = 5 % 



Channel 


A Hm (TeV) 

(va = - 1) 


A« m (TeV) 
iVij = + i) 


LL 


24.0 


16.4 


RR 


24.0 


16.5 


RL 


21.4 


17.6 


LR 


21.7 


17.4 
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TABLE VII. Expected 95% CL lower limits, A Jim , on 
the compositeness scale for different helicity channels of the 
quark-electron currents for L — 100 at 14 TeV with 

8e k = 15 % and 5L = 5 % 



Channel 


A iim (TeV) 

(vn = - 1) 


Aiim (TeV) 
(va = + 1) 


LL 


33.8 


20.1 


RR 


33.7 


20.2 


RL 


29.2 


22.1 


LR 


29.7 


21.8 



TABLE VIII. A5 CT for five different integrated luminosities 
for m, = - 1 at y/s = 14 TeV with Se k = 15 % and 
<5L = 5 % 



Asa 



Channel 


10 fb-' 


50 fb~' 


100 fb- L 


200 fb- 1 


500 fb'' 


LL 


16.0 


20.6 


23.4 


26.2 


31.0 


RR 


16.0 


20.5 


23.3 


26.2 


30.8 


RL 


15.1 


18.6 


20.9 


23.2 


26.8 


LR 


15.1 


19.1 


21.1 


23.5 


27.1 



TABLE IX. Asct for five different integrated luminosities 
for m, = + 1 at y/s = 14 TeV with Se k = 15 % and 
8L = 5 % 



A 5 q (TeV) 



Channel 


10 fb-^ 


50 fb-' 


100 fb-' 


200 fb-' 


500 fb-' 


LL 


12.4 


14.9 


15.9 


17.1 


18.3 


RR 


12.4 


14.9 


16.0 


17.1 


18.4 


RL 


13.1 


15.8 


17.3 


18.5 


20.2 


LR 


13.0 


15.7 


17.2 


18.3 


20.0 
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FIG. 6. 5<t discovery limit versus the integrated luminosity 
for rjij = — 1 (constructive interference). 
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FIG. 7. 5a discovery limit versus the integrated luminosity 
for rjij = + 1 (destructive interference). 
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